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Abstract Taking into account the helium and metal diffusion, we 
explore the possible evolutionary status and perform seismic anal- 
ysis of MOST target: the star e Eridani. We adopt the different 
input parameters to construct the models by fitting the available 
observational constraints: e.g., T e ff, L, R, [Fe/H]. From com- 
putation, we obtain the average large spacings of e Eridani about 
194 ± 1/iHz. The age of the diffused models has been found to 
be about 1 Gyr, which is younger than one determined previously 
by models without diffusion. We found that the effect of pure he- 
lium diffusion on the internal structure of the young low-mass star 
is slight, but the metal diffusion influence is obvious. The metal 
diffusion leads the models to have much higher temperature in the 
radiation interior, correspondingly the higher sound speed in the 
interior of the model, thereby the larger frequency and spacings. 

Key words: stars: evolution — stars: interiors — stars: individual: 
e Eridani 



1 INTRODUCTION 

Element diffusion, sometimes named "atomic" diffusion, is a basic physical el- 
ement transport mechanism which is driven by pressure gradients (or gravity), 
temperature gradients and composition gradients. The effects of helium diffu- 
sion on the structure of solar models were first studied by Noerdlinger dl9771 l. 
Later helioseismology provided an effective method to study this physical pro- 
cess (Guenther et al. [19931 Guen ther [T9941 Bahcall & Pinsonneault |1992bl 
Christensen-Dalsgaard et al. 1993; Guenther & Demarque |1997l l. Recently, he- 
lioseismic inferences have demonstrated that models incorporating element dif- 
fusion are in substantially better agreement with the inferred sound speed than 
the models that neglect it. 
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Up to date, the effect of element diffusion has been tested by asteroseismic 
method in solar type stars. According to the "second differences", Vauclair & 
Theado d2004b . Theado et al. (120051) and Castro & Vauclair d2006l presented 
discussion of asteroseismic signatures of pure helium diffusion in main-sequence 
stars between \.\Mq and 2.0M Q . In the present work, we study the effects of 
both helium and metal element diffusion on the models with mass less than 1M Q . 

Metal element could affect temperature gradient and convection in the enve- 
lope of stars. It also leads to the variations of the central temperature, density, 
pressure and radiative opacity. It is important to consider the influence of metal 
element on the internal structure of star. In this paper we mainly study the effects 
of metal diffusion on the stellar structure, evolution and the oscillation frequency 
for the late K-type star e Eridani. 

The bright K2V dwarf e Eridani (HD22049, HIP16537, HR1084, PLX742) 
is one of the nearest solar-like stars, with a distance of about 3.218pc. It has a 
planetary companion e Eridani b (Hatzes et al. |2000). Its proximity makes this 
planetary system as a prime target for future extrasolar planet direct - imaging 
efforts and studying the process of planetary systems formation. The success of 
these efforts will depend on the mass and the age of the system. Therefore the 
study of exoplanet - host star e Eridani is a key issue. 

In order to determine precise global parameters of the star e Eridani, some 
authors performed numerous theoretical analysis (Guenther & Demarque ll986l 
Guenther [19871 Soderblom & Dappen [19891 ( hereafter GD86, G87 and SD89 
respectively). Thanked for Noyes et al. (|19841 l, they measured the p-mode os- 
cillation spectrum of e Eridani and identified three individual frequencies in the 
power spectrum between 1500 fiHz and 2000 /iHz with average large spacing of 
172 ± 5 /iHz. Based on non-asteroseismic and asteroseismic observational con- 
strains, GD86, G87 and SD89 constructed a series of stellar models of e Eridani. 
However, in their works for e Eridani, there is a contradiction between the young 
age implied by the high chromospheric activity, rapid rotation rate, and the old 
age implied by the slightly metal-poor composition. In their works, GD86, G87 
and SD89 discussed this contradiction and given the age 12 Gyr, 10 Gyr and 1 
Gyr respectively. 

Recently many evidences from observations seem to indicate that e Eridani 
is a young main-sequence star with age less than 1 Gyr. For example, using Li 
abundances with the star's position in the H-R diagram and kinematics, Song et 
al. d2000) derived the age of 0.73 ± 0.2 Gyr; Measuring the radius of e Eridani 
by long-baseline interferometry, Di Folco et al. ([2004) estimated the age at 0.85 
Gyr; Saffe et al. (2005) used the calibrations of Donahue dl993l l and Rocha-Pinto 
& Maciel (1998) (which corrected the age with an effect from chromospheric 
activity) to estimate the ages of 0.66 Gyr and 0.82 Gyr respectively. SIMBAD 
database describes the age of e Eridani to be 0.66 Gyr. 

The aim of present paper is to provide the basic theoretical studies of e 
Eridani. We will use the latest interferometric observational radius (Di Folco 
et al. 2004) to constraint the stellar models and predict the frequency spacings. 
Then, we will investigate the effects of helium and metal diffusion on the internal 
structure, and the frequency spacings which also depend on the internal structure 
of the star. 

In section 2, the global parameters of the star were summarized. The input 
physics, computing method and results analysis are summarized in section 3. The 
helium and metal diffusion effects on the stellar models, and the oscillation fre- 
quency are discussed in section 4. Finally, the conclusion is presented in section 
5. 
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Table 1 The observational fundamental parameters of e Eridani. 



parameters 


e Hridani 


Ref 


M/Mq 


0.85 ± 0.04 


(1) 




0.80 


(2) 




0.82 


(3) 




0.83 


(4) 


[F e /H] sur f 


-0.13 ±0.04 


(5) 


R/ R Q 


0.743 ± 0.010 


(6) 


M v 


6.18 ±0.11 


(7) 




6.18 


(3) 


log(L/L ) 


-0.4616 ±0.044 


(8) 


T eff (K) 


5012 ± 67 


(9) 



(l)Santos et al. (20011; (2)Minier & Lineweaver (20061; (3)Takeda et al. ( [20051 ; 
(4) Fische r & Valenti J20051 : (5)Santos et al. ( 200 41; (6)D i Folco et al. ( [20041 ; (7)Song et 
al. J2000T >; (8)this paper; (9)Rami'rez & Melendez (2004) 



2 OBSERVATIONAL CONSTRAINTS 

Observationally it is difficult to determine accurately the mass of star except 
some binary systems, e Eridani is a single star and it has various data of masses 
which were published in the literatures, for example 0.85 ± O.O4A/ (Santos et 
al. p2"00"TT >, O.8OM (Minier & Lin eweaver [2006b , O.82M (Takeda et al. 120051 , 
O.83M (Fischer & Valenti [2005>. 

For the effective temperature of e Eridani, we adopt T e f / = 5012 ± 67 which 
was determined from the infrared flux method (IRFM) (Ramirez & Melendez 
120041 

The luminosity can be calculated using the absolute visual magnitude M v = 
6.18 ±0.11 (Song et al. [2000l Takeda et al. 120051 . Given the bolometric correc- 
tion B.C. = - 0.28 (Takeda et al. 2005), the absolute bolometric magnitude is 
therefore Mt, i = M v + B.C. = 5.9 ±0.11. Then the luminosity is computed by 
formulae: Mboi = - 2.5log(L / X ) + M& o i i0 . Adopting Mboi,o = 4.746, we thus 
obtain log (L/L Q ) = - 0.4176 ~- 0.5056'. 

e Eridani is slightly metal-poor compared with the Sun. Santos et al. (2004) 
obtained [Fe/H] = - 0.13 ± 0.04 from the spectroscopic analysis basing on 39 
Fei and 12 Fen lines. In order to deduce [Z/X], we assume that the logarith- 
mic value of \Zj X] is proportional to the abundance ratio [Fe /H] (Thoul et al. 
120031 : 

log[Z/X] star = [Fe/H] star + log[Z/X] Q (1) 

where [Z/X]@ = 0.023 (Grevesse & Sauval [T998l l. Using this ratio, we obtain 
[Z/X] surf = 0.017 ±0.0016. 

A star's initial helium mass fraction Yi is important for determining its struc- 
ture, but Yi of e Eridani is unknown by observations. GD86, G87 adopted Yi to be 
0.236 and SD89 considered Yi to be from 0.24 to 0.26. In our work, we reference 
the value given by GD86, G87, SD89 and adopt the initial helium abundance Yi 
= 0.245 ±0.015. 

The frequency, especially large spacings, sensitively depend on the radius. In 
order to predict the large spacings which will be observed by MOST, we take the 
latest interferometric observational radius R/Rq = 0.734 ± 0.01 (Di Folco et al. 
120041 . 
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3 STELLAR MODELING 

3.1 Input physics 

The evolutionary models are computed using the Yale stellar evolution code ver- 
sion (Guenther et al. 1 19941 1 which has been modified to include the effects of 
element diffusion. The initial zero-age main sequence (ZAMS) models were cal- 
culated from pre-main sequence evolution which were assumed to have spheri- 
cal symmetry, no rotation, no magnetic field. In the computation we use OPAL 
equation of state tables EOS2001 (Rogers & Nayfonov 2002 ), the opacities with 
a smooth blend of OPAL GN93 (Iglesias & Rogers [19961) and OPAL tables 
(Alexander & Ferguson 1994 ). The relevant nuclear reaction rates and cross sec- 
tions are from Bahcall & Pinsonneault dl992al |1992bl ) and Bahcall (I19891 ). e 
Eridani is a solar-like star, we choose the Krishna-Swamy ( 1 19661 ) atmosphere 
model. The diffusion of both helium and metal element abundance was consid- 
ered in the stellar model computation, by using the coefficients of Thoul et al. 
( 119941 ). In the next section we will introduce it in detail. 

Using the standard mixing-length theory, we set a = 1 .7 for all models, close 
to the value which is required to reproduce the solar radius under the same 
physical assumptions and stellar evolution code (Murphy & Demarque 2004). 
If rotation and magnetic field are neglected in the construction of stellar model, 
convection is the most important mechanism that is usually included in stan- 
dard stellar models which describes material motions within a star especially for 
low-mass main sequence star. The convective overshoot is an another factor to 
influence the chemical composition of the star. The change of mixing-length pa- 
rameter a would also influence principally the radius of the base of the external 
convection zone. The presence of core overshoot would extend the core-burning 
phase of evolution and increase the estimated age of star. This important topic 
has been studied separately by Bi et al. (2008 ). However, the element diffusion 
could also lead to the redistribution of the element in star and change the depth of 
the convective envelope (Guenther 1994). In order to test the effect of diffusion 
on the chemical composition, internal structure and the frequency of e Eridani, 
we would not discuss the influence of the convection and overshoot in present 
work. 

3.2 Element diffusion 

The element diffusion in stars is driven by gravitational setting, temperature gra- 
dients, composition gradients and radiation pressure. It is described by Bahcall 
& Pinsonneault ( 1992a ) in detail. The gravity and temperature gradients tend to 
concentrate the helium and heavier elements toward the center of the star, while 
the hydrogen diffuses outward. The concentration gradients oppose the above 
processes. In addition, radiation pressure can cause partially ionized or neutral 
species to rise relative to species with a small cross section. However, the radi- 
ation pressure are efficient in the external regions of main-sequence stars with 
T e ff > 6000i4T (Michaud ll9761 ). So the radiative pressure causes negligible dif- 
fusion in the K-type star and will be neglected in this paper. 

3.2.1 Basic equation for element diffusion 

Most works on chemical diffusion use either the Chapman - Enskog procedure 
(Chapman & Cowling 119701 ) or the method of Burgers (1969 ) for deriving the 
transport properties from the Boltzmann equation. In our work the element dif- 
fusion is described using the equations of Burgers (119691 ) with diffusion velocity 
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coefficients from Thoul et al (1994). A single model includes both helium and 
metal element diffusion. The change rate of the element mass fractions due to 
the diffusion is now written as (Thoul et al l 19941 Bahcall & Loeb ll990l l: 

dX s = 1 d[r*X s T 5 /%(r)} 
dt pr 2 dr 

where the partial derivatives are evaluated in the local rest frame of a mass shell 
in the star, p, r, T are the local value of density, radius, temperature and X s is 
the mass fraction of element s. 

The diffusion velocity for species s is defined by 



(3) 



and the function £ s (r) is expressed as 



In equation (4), A p , At, A c correspond to gravitational setting, thermal diffusion 
and concentration gradient diffusion which are the functions of the mass fractions 
only. Species c can equal to s, but c aren't electron and helium element. 

3.2.2 The treatment of diffusion in the evolutionary code 

The Yale stellar evolution code has been modified to include the effects of 
He diffusion from gravitational settling and thermal diffusion (Bahcall & 
Pinsonneault 1992a) using the method of Bahcall & Loeb ( 1990 ) and added the 
metal diffusion in 1994. About the method by which treatment of element dif- 
fusion in a standard stellar evolution code, it has been described by Bahcall & 
Pinsonneault (1992a) in detail. In this section, we summarize it again. 

In the diffusion subroutine written by Bahcall and Pinsonneault, the main 
problem is to solve the diffusion equation and calculate the change of the element 
abundance by diffusion. The diffusion equation (2) can be written as following 
form: 

dx 1 . d . ^ d dx 
dt pr 2 dr dr dr 

where D\, D2 is the diffusion coefficients. 
For helium diffusion 



2 ^.W + -(D 2 — )] (5) 



Di(Y) = Fay 'CT' 2 -^-X-iAX+AX) 



For metal diffusion 




r,5/2 



d In p 
In A _ ,_ dr 
Z 



F ^. r 2. T 5/2 

In A 



A 



(6) 



(7) 



and radius, temperature, and density are in nondimensional units defined by r = 
r'/R Q , T = T'/10 7 K, p = p'/l OOg cm " 3 and t = t'/10 13 yr (Bahcall & 
Pinsonneault |1992al Chaboyer et al. [I^9"2"i ). 

In expressions (6) and (7), F gy and F gz are the adjustment factors by which 
the helium and metal diffusion coefficients are multiplied. In our work, we set 
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both of the two factors to 1.0. A p , At, A c are the diffusion velocity coefficients 
as described in equation (4). The exact numerical solution of these coefficients 
are calculated by an export subroutine which was developed by Thoul in 1994. 
The computation method is exhibited detailedly by Thoul et al. ( 1994). 

In equation (5) the partial time derivative is evaluated at constant mass shell 
of the star. The diffusion equation is solved with zero hydrodynamic velocity of 
the stellar plasma. The right hand of diffusion equation (5) has two terms: the 
first term depends on first spatial derivatives of the element mass fraction and 
the second term depends on second spatial derivatives of the element abundance. 
The first term is solved explicitly using the two-step Lax-Wendroff technique 
(Press et al. 1986 ) by neglecting the second derivatives of the element abundance. 
We then use this trial solution as the initial abundance to determine the second 
derivatives by a fully implicit method (Bahcall & Pinsonneault 1992a). At last, 
the element abundance were updated. 

For the helium diffusion, we only consider 4 He diffusion and ignore 3 He 
diffusion (Loeb et al. |1989l ). The time rate of change of the 4 He mass fraction Y 
is equal in magnitude and opposite in sign to the rate of change of the hydrogen 
mass fraction. 

The diffusion of all metal elements were assumed to diffuse at the same rate 
as fully ionized iron. Because of the diffusion of metal elements, the radiative 
opacity have to be changed at each spherical shell in the stellar model after each 
time step. We calculate the effect of metal element diffusion on the opacity by 
computing a total metal element abundance at each model radius Z(r), and then 
interpolating for the opacity between opacity tables with different total metal 
element abundances (Bahcall et al. 119951 ). 

The diffusion subroutine carries out the diffusion calculations using data sup- 
plied by other parts of the Yale code, in which the thermal structure and the 
element abundances are calculated. It is assumed that the amount of diffusion 
within a given time step was too small to affect significantly the changes in ther- 
mal structure, in abundances, and in nuclear reaction rates that are calculated 
elsewhere. After the diffusion computation, we update the element abundance 
and then calculate the thermal structure. 

In the diffusion subroutine, the element diffusion is treated only in the radi- 
ation region. The characteristic time for element to diffuse a solar radius under 
solar conditions is of the order of 10 13 yr (Bahcall et al. 1 19951 ). much larger than 
the age of the sun and e Eridani. However, helioseismic inferences have demon- 
strated the significance of the element diffusion (Cox et al. 1989 1 Bahcall & Loeb 
119901 Bahcall & Pinsonneault |1992al Proffitt [T9"9"4"i ), we would test how much ef- 
fect on the internal structure, the age and the frequency of e Eridani by both the 
helium and metal diffusion using the asteroseismic method. 

3.3 Computational method and results analysis 

3.3.1 Evolutionary tracks and candidates for pulsation 

In order to reproduce observational constraints of e Eridani, we have computed 
a grid of evolutionary tracks for six masses M/M Q : 0.80, 0.81, 0.82, 0.83, 0.84, 
0.85 and the initial metallicity Z t : 0.011, 0.012, 0.013, 0.014. 

From a series of evolutionary tracks which are computed, we select those 
models that land within the observational error box (being composed by 
logT e //, logL/L©, log i?/i?©) in the theoretical H-R diagram. Then we fur- 
ther choose the models on the basis of the constraint of the observational [Z/X]. 
In terms of the computation alone, we note that the age of e Eridani is likely in 
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Table 2 Model Parameters of e Eridani. 



Model 


Al 


A2 


A3 


Bl 


B2 


B3 


CI 


(J2 


Ci 






M/Mq = 


0.83 




M/Mq = 


0.83 




M/Mq = 


0.83 






Zi = 0.012 




Zi = 0.013 




Zi = 0.014 






Yi = 0.2315 




Yi = 0.2365 




Y t = 0.2415 


Dif. 


None. 


Y.Dif 


Y& Z.Dif 


None. 


Y.Dif 


Y& Z.Dif 


None. 


Y.Dif 


Y& Z.Dif 


Age(Gyr) 


1.151 


1.089 


1.0 


1.206 


1.138 


1.0 


1.264 


1.191 


1.0 


L/L e 


0.321 


0.321 


0.321 


0.319 


0.319 


0.319 


0.317 


0.317 


0.317 


R/Rq 


0.749 


0.750 


0.744 


0.750 


0.751 


0.743 


0.751 


0.753 


0.741 


T eff 


5024.9 


5019.8 


5039.5 


5012.2 


5006.7 


5036.3 


4999.9 


4993.8 


5033.5 


Tcz/R 


0.7073 


0.7053 


0.7048 


0.7046 


0.7033 


0.7024 


0.7024 


0.7022 


0.7000 


Tcz(s) 


1504.8 


1514.1 


1496.7 


1514.4 


1522.9 


1496.7 


1523.2 


1529.5 


1497.0 


T (S) 


2524.5 


2532.5 


2502.4 


2530.1 


2538.8 


2493.6 


2535.7 


2545.2 


2484.5 


Tcz(s)/t (s) 


0.5961 


0.5979 


0.5981 


0.5986 


0.5999 


0.6002 


0.6007 


0.6009 


0.6025 


< Av> 


192.2 


191.4 


193.8 


191.6 


191.0 


194.6 


191.2 


190.5 


195.1 



the range 0. 1 Gyr - 1 .0 Gyr. Here we consider the age (< 1 Gyr) as a property of 
star which be same as SD89. According to the age estimation from observation, 
which given in the introduction, we constraint the age of models in the range 
0.60 Gyr - 1.0 Gyr. At last, 12 tracks fall within the observational error box and 
are plotted in Fig. 1 (a). 

In order to deduce the set of parameters lending to the better agreement with 
observations, we perform a \ 2 minimization method described by Eggenberger 
(2005). We define the % 2 function as follows: 



X 2 s £( - f (8) 

i=l 



where the vector C° bs contains the observations: C° bs = (L/Lq, T e ff, R/Rq, 
\Z/H] S ). The value C° bs and the vector a which contains the errors on the ob- 
servations are given in Table 1. From this method we obtain three models A3, 
B3 and C3 which are given in Fig. 1 (b) and Table 2. 



3.3.2 Calibration of models with and without diffusion 

In order to investigate the effect of element diffusion, especially metal diffu- 
sion, basing on the selected models A3, B3 and C3, we calibrate three groups 
of models with mass M/Mq = 0.83. There are the different initial chemical 
compositions among the different groups. The models at the same group have 
the same initial chemical compositions, but they have different types of element 
diffusion. The detail description of the models are given in Table 2. 

We use the stellar pulsation code of Guenther (119941 ) to calculate the eigen- 
frequencies and the large spacings of models listed in Table 2. For solar-like 
stars, the eigenfrequencies u n j, of oscillation modes, characterized by the radial 
order n at harmonic degree /, satisfy the simplified asymptotic relation (Tassoul 
rT980b : 

v n .i = Ai/(n + - + a + -) + e n , z (9) 

The large spacings is defined by the frequencies of the same harmonic degree / 
and the adjacent radial order n: 

AL> n< l = v n> i - v n -x t i (10) 

We average the large spacings over the modes 1 = 0, 1, 2, 3 and n = 10, 11, 30 
to obtain the average large spacings. 
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-0.6 [ 1 

3.730 3.720 3.710 3.700 3.690 3.680 

Log(Teff) 





-0.40 


(b) 







J 


-0.45 


























Logl 


-0.50 






A3 

B3%X3 








-0.55 











3.708 3.704 3.700 3.696 3.692 

Log(Teff) 

Fig. 1 (a)Evolutionary tracks (12 in total) in the H-R diagram from 
ZAMS. The solid lines indicate the boxes delimited by the observed lu- 
minosity and effective temperature. The two dash dot lines denote the 
boxes delimited by the interferometric radius. (b)Enlarge the error box 
of Fig (a). The model A3, B3 and C3 are denoted by filled circle, triangle 
and five pointed star respectively. 



From the age of the star in Table 2, we can see that the diffusion could speed 
up the evolution of stars. The model with both helium and metal diffusion has the 
smallest age of all the models in the same group. The average large spacings is 
about 194 ± 1 /iHz of models A3, B3 and C3. The pure helium diffusion hardly 
affect v n _i, but the effect is relatively obvious after adding the metal diffusion. 
It also seems that the diffusion effect is increase alone with the increase of the 
initial heavy element abundance (Yi, Zi). In the following section, we select 
models in group (C) to discuss differences in their internal structures in detail 
using the asteroseismic method. 

4 ASTEROSEISMIC TEST HELIUM AND METAL DIFFUSION 
EFFECTS 

4.1 "Second differences" tests helium and metal gradients 

Stellar acoustic p-modes with low I degree can propagate deeply inside the stars. 
However, as mentioned by Gough (1990), rapid variations of the sound speed in- 
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(a) Zi = 0.014 
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2500 2000 1500 1000 500 
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2500 2000 1500 1000 500 

t(s) 
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0123456789 10 

t(1000s) 



Fig. 2 The stellar structures of the group (C) models in Table 2 with 
= 0.2415, Zi = 0.014. The detailed notes are labeled in the figure. Solid 
lines: standard homogeneous model CI without diffusion; dashed lines: 
model C2 with pure helium diffusion; dash dot lines: model C3 with 
both helium diffusion and metal diffusion. The panel (a) helium profiles 
as a function of the acoustic depths in the models, i.e. the time for the 
acoustic waves to travel from the surface down to the considered layer; 
(b) the metal profiles: a clear gradient at the bottom of the convection 
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side a star lead to partial reflections of the sound waves. A conveniently and eas- 
ily evaluated measure of this oscillatory component is the "second differences" 
with respect to the radial order n of the frequencies v n i for the same value of the 
harmonic degree I: 



This measure is contaminated less than the first difference Aiu n i = v n .i - v n -i t i 
by the smoothly varying components of Aii/ n ;. This modulation in the frequency 
has been extensively studied in the Hen ionization zone and the edge between 
the convective and radiative zones. In the Hen ionization zone there is a rapid 
variation of the adiabatic exponent Ti and in the base of the convection zone 
there is essentially a discontinuity in the sound speed gradient, i.e. dc/dr, pro- 
duce a discontinuity in v n \. In order to identify the different components which 
modulate the oscillations, we have computed the Fourier transform of the "sec- 
ond differences". The modulation period of the "second differences" is twice the 
"acoustic depth" of the region where the feature occurs: 



where r s is the time needed for the acoustic waves to travel between the surface 
and the considered region (the Hen ionization zone or the base of the convection 
zone), c(r) is the sound speed at radius r, and r s is the radius of the considered 
region. 

Vauclair & Theado (l2004t . Theado et al. (|2005T > and Castro & Vauclair (12006b 
described asteroseismic signatures of pure helium diffusion in stars between 
1 . 1 Mq and 2.0A/ Q in terms of "second differences". They have studied the pre- 
cise signatures on the oscillation frequencies of helium gradients inside stars. In 
our work we compute stars of O.83M in which not only including helium diffu- 
sion but also containing metal diffusion. The aim of the present section is to study 
the precise signatures on the oscillation frequencies of not only helium gradients 
but also metal gradients inside stars and to test the effect of metal diffusion on 
the internal structure. 

Fig. 2(a)(b) display the helium and metal abundance profiles in the mod- 
els without element diffusion (solid line), with pure helium diffusion (dashed 
line) and with both helium and metal diffusion (dash dot line) as a function of 
the acoustic depth. Due to diffusion, helium and metal drift inward and form a 
gradient just below the convective zone. The acoustic depth of the bottom of the 
convection zones t cz and the corresponding t cz /tq , where To is the total acoustic 
depth of the star, are listed in Table 2 about models with different initial chem- 
ical composition and different types of diffusion. Fig. 2(c)(d) give the adiabatic 
exponent Ti profiles and the sound speed gradients dc/dr versus the acoustic 
depth. From the dc/dr profiles it is clear to see that the dips around 300s and 
1500s are caused by the Hen ionization zones and the bottom of the convection 
zones. 

Considering the asymptotic approximation validity, we select modes of de- 
grees 1 = 0, 1, 2, 3 and oscillation frequencies between 1500/iHz and 3000fj,Hz 
to compute the "second differences". In Fig. 2(e)(f) we present the second dif- 
ferences of our models and the Fourier transforms of these curves. We have ex- 
plained that the modulation period corresponding to the peaks in the Fourier 
transforms are twice the "acoustic depth" of the region where the feature occurs. 
By comparing the dips of dc/dr profiles in Fig. 2(d) with Fourier transforms in 
Fig. 2(f), one can easily recognize the peaks due to the Hen ionization zones and 
those due to the bottom of the convective zones in Fig. 2(f). 



AlVni = V n -l,l — + V n +l,l 



(ID 




(12) 
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From the Fourier transforms in Fig. 2(f), we can see that the amplitudes of 
the peaks corresponding to the modulation periods due to the base of the convec- 
tion zones of models CI, C2, C3 are different. Mazumdar & Antia ( 1200 11 1 shown 
that the amplitude of the oscillatory signal in the second differences contains an 
amplification factor of 4 sin 2 (7rr s /ro). When t s /to > i, the Fourier transforms 
show the peak amplitudes decrease for deeper layers (Vauclair & Theado 2004). 
In terms of r cz /R and t cz /tq in Table 2, we find that diffusion could deepen 
the outer convective zone, especially including the metal diffusion, the base of 
the convective zone is deepest of all in the same group. From the Fourier trans- 
forms, we can see that the peak amplitude at 2t cz of models with pure helium 
diffusion is lower than that without diffusion. In term of amplification factor of 
4 sin 2 (7TT s /ro), the peak amplitude of models with both helium and metal diffu- 
sion should be lowest of all, but it is higher than that with pure helium diffusion. 
So we think that the peak amplitude is mainly determined by the depth of the 
base of the convection zone, meanwhile it is also sensitive to the element gradi- 
ent at the bottom of convective zone. From Fig. 2(b), it is clear to see that the 
metal gradient steepen more rapidly than helium gradient. There is a stronger 
reflection of the sound waves in the region of the metal gradient, that is why the 
model contain metal diffusion have a higher amplitude of the peak of r cz . 



4.2 Tests of internal structure 



In this section, we compute the so-called small frequency spacings and frequency 
separation ratios of the low-/ p-modes to test the internal structure of models in 
which considered different types of diffusion. The small frequency spacings is 
given by the combination: 

du+2(n) = v n ,i - v n -i : i+2 (13) 

Using the asymptotic theory of p modes it can be shown that (Christensen- 
Dalsgaard & Berthomieu fTWTl Basu et al. l2007l : 

du+2{ n)c-(Al + 6)^L (14) 
47T - z v n ,i Jo ar r 

where R is the stellar radius and Ai/„; is the large spacings. In the core the 
gradient of the sound speed is large and r is small, the integral in equation (14) is 
dominated by conditions in the core. So the small spacings usually test the stellar 
core. But the small spacings are also slightly affected by near-surface effects. In 
order to reduce the effect of near surface uncertainties, the frequency separation 
ratios are used (Roxburgh & Vorontsov 2003): 

/ \ d 02 (n) d 13 (n) 
Ai(ra) A (n + 1) 

Each panel of Fig. 3 show small spacings rfo2 and frequency separation ratios 
ro2 for different models in group (C) with different types of diffusion which are 
labeled in the figure. From the small spacings and frequency separation ratios we 
can see that the pure helium diffusion hardly alter the internal structure in these 
young low-mass models. By comparing the models included metal diffusion 
with the models no element diffusion, the differences of their internal structures 
are relatively obvious. 

In order to explain the changes of the internal structure from the astero- 
seismic test, we compare the internal physical parameters of models with and 
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Fig. 3 Upper part: small spacings do2(n), lower part: frequency sepa- 
ration ratios r 02 (n) for the models in group (C). The detailed notes are 
labeled in the figure. 



without diffusion. It is known that the sound speed depends on both the mean 
molecular weight and the temperature (Basu et al. 2007; Bi et al. 120081 ) 

2 KbT T 
c ~ oc — (16) 

In Fig. 4, we give the 8c/ c, 5T/T, and S/i/ fi variations as function of radius. 

From Fig. 4, it is clear to see that the effect of the pure helium diffusion on 
the sound speed c, temperature T and the mean molecular weight fi is small. The 
helium diffusion hardly alter the c, T and /i in the radiation interior and the effect 
mainly concentrate on the convection zone. Due to diffusion, helium drift inward 
just below the convection zone. It induce /i is smaller than that without diffusion 
in the envelope. On the other hand, the increase of the opacity in the envelope due 
to the settling of He leads to the lower temperature than that without diffusion. 
However, the helium diffusion decrease both the T and fi, the sound speed c 
slightly in the convective envelope. So the small spacings and separation ratios 
are similar with and without helium diffusion. 

In Fig. 4 the dash dot lines describe the differences of c, T and /i between the 
model with both helium and metal diffusion and the model without any diffusion. 
From Fig. 4, we can see that the effect of metal diffusion on the c, T and /i not 
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only in the convective envelope but also in the radiation interior. Due to diffusion, 
helium and metal drift inward just below the convection zone and decrease fi in 
the outer envelope. At the same time, the depletion of Z increase the opacity in 
the envelope, otherwise the concentration decrease the opacity in the radiation 
interior. It leads to lower temperature in the envelope and higher value in the 
interior. From Fig. 4, we can see that the temperature is more sensitive to the 
metal diffusion than other physical parameters. The Sc/c is mainly dependent 
on the variation of 5T/T. Except the convective envelope, the sound speed c is 
much larger in the interior of the model. From the expression of the sound speed 
gradient in terms of the equation (16): 

dc = _i r ^ l^dT 

dr dr 2 V dr 2 (Tfi)i dr V ' 

we find that the temperature T is the dominate factor in the sound speed gradient 
in causing the differences between the models with and without metal diffusion. 
Due to dT/dr is negative, \dc/dr\ is larger than that of model without metal 
diffusion in the radiation interior. So from equation (14), the small spacings are 
larger than other models. That are explained the variation of the internal structure 
including metal diffusion. 



5 CONCLUSIONS 



Taking into account the effects of helium and metal diffusion on stellar evolution, 
we have presented a detailed modeling of the MOST target: the star e Eridani by 
fitting the available observational constraints: effective temperature, luminosity, 
high precision interferometric observational radius and metallicity. We also con- 
sider the age (< 1 Gyr) as a property of star to constraint the models, which are 
the same as SD89. 

According to our results, the location in the H-R diagram indicates that the 
star e Eridani has evolved little since its arrival on the zero-age main sequence. 
Taking into account the helium and metal diffusion, we find that the mean large 
spacings < Av >, which are averaged over 1 = 0, 1, 2, 3 and n = 10, 11, 12, 
30, is about 194 ± 1 /j,Hz. 

In order to test the effect of helium and metal diffusion, we construct three 
groups of models with the same mass but different Yj and Zi. Then, for the 
given chemical composition we calibrate three models with the different types 
of diffusion which are listed in Table 2. In our work we fixed the mixing length 
parameter a = 1 .70 and neglected convective overshoot. The results show that 
the age of diffused models is about 1 Gyr, and be younger than that determined 
by the models without diffusion. 

Due to the effect of diffusion, helium and metal fall below the outer convec- 
tive zone and form a gradient. We use the "second differences" to test the depth 
of the convective zones. When metal diffusion is considered, the convective zone 
deepens more quickly and the metal gradient steepens more rapidly than the he- 
lium gradient. That induce a stronger reflection of the sound waves just below 
the convection zone, which explains the higher amplitude of the peak related to 
the base of the convective zone than that of pure helium diffusion. 

The small spacings and frequency separation ratios have been used to test the 
differences of internal structure in models with and without diffusion. According 
to our results, the pure helium diffusion hardly alter the internal structure in 
the young low-mass models. Taking into account the metal diffusion, it mainly 
induce a much higher temperature in the radiation interior, thereby the higher 
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sound speed in the interior of the model. Correspondingly, the frequencies and 
large spacings are larger obviously than that of models without metal diffusion. 
The variation of internal structures between models with and without metal dif- 
fusion is increase obviously alone with the increase of initial metallicity. So we 
conclude that it is necessary to consider the metal diffusion in the young main- 
sequence low-mass star. 
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Fig. 4 From up to down, each panel of figure give the relative sound- 
speed difference 5c/ c, the temperature difference 8T/T and the mean 
molecular weight <5^//i respectively for the models in group (C). All 
differences are with respect to the model CI which without any diffusion. 
The dashed lines render the differences between the model with pure 
He diffusion and the model without any diffusion (e.g. 5fi/fi= (/J-y.Dif 
- UNone.Dif) I /^None.Dif)- The dash dot lines render the differences 
between the model with He & Z diffusion and the model without any 

diffusion (e.g. 8[lj\l = (py&Z.Dif - ^None.Dif) I ^None.Dif)- 



